Interplay of complete wetting, critical adsorption, and capillary condensation. 
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The excess adsorption F in two-dimensional Ising strips (oo x L) subject to identical boundary 
fields, at both one-dimensional surfaces decaying in the orthogonal direction j as —hij~ p , is studied 
for various values of p and along various thermodynamic paths below the critical point by means 
of the density-matrix renormalization-group method. The crossover behavior between the complete 
wetting and critical adsorption regimes, occurring in semi-infinite systems, are strongly influenced 
by confinement effects. Along isotherms T = const the asymptotic power law dependences on the 
external bulk field, which characterize these two regimes, are undercut by capillary condensation. 
Along the pseudo first-order phase coexistence line of the strips, which varies with temperature, we 
find a broad crossover regime where both the thickness of the wetting film and F increase as function 
of the reduced temperature r but do not follow any power law. Above the wetting temperature the 
order parameter profiles are not slab-like but exhibit wide interfacial variations and pronounced 
tails. Inter alia, our explicit calculations demonstrate that, contrary to opposite claims by Kroll 
and Lipowsky [Phys. Rev. B 28, 5273 (1983)], for p — 2 critical wetting transitions do exist and we 
determine the corresponding wetting phase diagram in the (hi,T) plane. 

PACS numbers: 05.50.+q, 68.35.Rh, 68.08.Bc 



I. INTRODUCTION 



Far away from phase boundaries, i.e., deep in the one- 
phase region, condensed matter is perturbed by confining 
walls only within a thin layer proportional to the bulk 
correlation length £, i.e., at most approximately 15 A. 
This changes drastically if the thermodynamic state of 
the bulk system is moved towards the boundary of phase 
transitions between the bulk phases. If the phase tran- 
sition is of second order the bulk correlation length di- 
verges and so-called critical adsorption occurs, i.e., the 
perturbation due to the wall located at y = pen- 
etrates deeply into the bulk, resulting in an algebraic 
divergence of the thickness of the interfacial structure 
Jo" (P(y) ~ Pb) dy/( P (0 + ) - p b ) ee T/(p(0+) - Pb ) ^ oo, 
where for a one-component fluid as a paradigmatic case 
Pb is the bulk number density for a given temperature 
T and chemical potential p, and p(y) is the fluid num- 
ber density profile [j], Q. Moreover, the order-parameter 
profile m(y) = (p(y) — p c )/Pc where p c is the bulk crit- 
ical number density, is governed by a universal scaling 
function which decays exponentially on the length scale 
of £ in the direction y normal to the surface [l], If 
the phase transition is of first order wetting phenom- 
ena occur as a result of a subtle interplay between the 
substrate potential, interactions among fluid particles, 
and entropic contributions, in particular interfacial fluc- 
tuations [3(. The thickness of the interfacial structure 
diverges too, but here with exponents which are deter- 
mined by the decay exponents of the pair potentials be- 
tween the fluid particles and of the substrate potential [|[ . 



These two types of phenomena are of significant practical 
importance, in particular for the condensed phase being 
a fluidum near the gas-liquid transition, such as water 
near the thermodynamic coexistence with its vapor, or 
binary liquid mixtures such as water and hydrocarbons 
near liquid-liquid phase transitions. In this context var- 
ious applications arise ranging from the use of colloidal 
suspensions Q to petroleum recovery Wetting films 
are relevant in many types of liquid coating processes, 
such as lubrication and adhesion and also for microflu- 
idics and nanoprinting @, 0] • Critical adsorption plays 
an important role, e.g;, for heterogeneous nucleation in 
supercritical solvents [8| or in micro- and nanofluidic sys- 
tems in order to achieve wetting of these small structures 

a- 

A characteristic property of liquids is that first-order 
gas-liquid or liquid-liquid phase transitions end at critical 
points. Therefore critical adsorption and wetting must 
emerge from each other upon moving along the first-order 
phase boundary. This creates a puzzle. Liquids governed 



by dispersion forces (decaying asymptotically < 



with 



distance r) belong to the Ising universality class [Tcj, i.e., 
the critical exponents describing the singular behavior of 
various thermodynamic quantities and structural prop- 
erties are those of the Ising universality class. This also 
holds at interfaces. Specifically, along bulk coexistence 
Am = (/i — /.io(T))/(fcsT c ) = one has asymptotically 



A/i = 0, 



(1) 



where r = (T — T c )/T c is the reduced deviation from 
the critical temperature T c and A/i measures, in units 
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gas 



FIG. 1: Schematic drawing of the liquid-gas coexistence curve 
A/i = in the ( Afi, T) bulk phase diagram. The thick solid 
line indicates the bulk coexistence line. T w is the transition 
temperature for a critical wetting transition. Thermodynamic 
paths (I) and (II), along which the behavior of the adsorption 
is discussed in the main text, are displayed as well as various 
crossover lines. The hatched area is the complete wetting 
regime. Note that 2>v ~ 1.89 is larger than A ~ 1.56. 



of ksT c , the deviation of the chemical potential from its 
value no(T) at two-phase coexistence, whereas along the 
path t = 



r~|A/*l 



(/3-^/A 



r = 0, 



(2) 



where the critical exponents v, /?, and A take the val- 
ues of the Ising model (in d = 3: v = 0.6301(4),/? = 

0. 32653(10), and A = 3u - (3 = 1.564 p]]). The lat- 
ter universal singularity T ~ |A/i| -0194 is weaker than 
the non-universal one V ~ |A/i| -1 / 3 for complete wetting 
3, although according to the general renormalization 
group (RG) arguments the dispersion forces, which are 
responsible for this non-universal behavior, at first glance 
should give rise only to corrections to scaling (i.e., sub- 
dominant power laws). Naturally the question arises how 
the system manages to restore universality (i.e., dom- 
inance by the weaker universal power law) upon mov- 
ing from the non-universal complete wetting behavior to 
the universal critical adsorption behavior. This issue has 
been addressed in Refs. [IH, [HI, [l4[ and the proposed 
scenario is such that upon approaching bulk coexistence, 

1. e., for t fixed and A/j, — ► 0~ (see path (II) in Fig.[l}, the 
system will always cross over ultimately to the complete 
wetting regime but such that upon approaching T c (i.e., 
for smaller |r|) this crossover occurs closer and closer to 
the coexistence curve. The proposed mechanism is based 



on the argument that long-ranged dispersion forces are 
relevant and can dominate if the thickness I ~ |A/.i| -1 / 3 
of the wetting layer is much larger than the bulk cor- 
relation length £ ~ |t| ~ v . As described more closely 
in Sec. [TTJ for dispersion forces in d = 3, this occurs if 
|A/z| <C |t| 3i/ (see the hatched area in Fig. [T|). Accord- 
ingly, upon approaching T c this complete wetting regime 
becomes vanishingly small. This expectation follows sim- 
ilar arguments as how dispersion forces become irrelevant 
for, e.g., two-point correlation functions in the bulk or for 
the critical adsorption profile itself, for which RG theory 
predicts exponential decays which, however, are domi- 
nated asymptotically by the algebraic decay due to the 
dispersion forces. Nonetheless the universal scaling func- 
tions, which capture the exponential decay, are expected 
to describe also the neighborhood of T c correctly, even in 
the presence of the irrelevant dispersion forces. 

So far this puzzle has not been resolved satisfactorily 
by explicit calculations. The main reason is that mean 
field theory (MFT) as the natural theoretical starting 
point does not capture the correct d = 3 exponents, 
which is essential for disentangling the different contri- 
butions and to describe the crossover. Within MFT, 
v = j3 and the adsorption diverges only ~ In |A/x|. In sys- 
tems governed by short-ranged forces only, within MFT 
both critical adsorption and complete wetting follow a 
divergence <~ In | A/i| and therefore they cannot be distin- 
guished; but the divergences differ beyond MFT. How- 
ever, going beyond MFT is a difficult theoretical chal- 
lenge, due to the spatial inhomogeneity. 

There is also an urgent experimental need to resolve 
this issue. An ellipsomctry study of gravity-thinned com- 
plete wetting layers in the binary-liquid mixture of cyclo- 
hexane and methanol [l5[ reports results which are not 
in accordance with the theoretical predictions. The data 
indicate the divergence of the film thickness upon ap- 
proaching the bulk critical temperature characterized by 
a critical exponent which is distinct from the expected, 
critical exponent for the bulk correlation length. Scaling 
arguments have been put forward according to which the 
observed effective exponent is associated with a broad, in- 
termediate scaling regime facilitating the crossover from 
complete wetting in the presence of dispersion forces to 
critical scaling [l5(. Also recent neutron reflcctometry 
data for the adsorption from alkane-pcrfluoroalkane mix- 
tures at fluorophobic and fluorophilic surfaces [l6| are not 
in agreement with the theoretical predictions. These au- 
thors have found that the behavior of the adsorption, as 
a function of temperature in the one-phase region upon 
approaching liquid-liquid coexistence with the composi- 
tion well removed from the critical composition, can be 
represented by a power law with an exponent which dif- 
fers from both that for complete wetting and that for 
critical adsorption [l6| . 

The predictions described above for the crossover be- 
tween the complete wetting regime of fluids with long- 
ranged interactions and the critical regime are based on 
general scaling arguments and on the analysis of simple 
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models for the effective interface Hamiltonian which are 
supposed to describe the relevant physics at length scales 
much larger than the bulk correlation length [13j . Only 
few micro scop ic studies are available which test these 
ideas 0, Gj . Those are based on density functional 
methods and are of mean-field character, which in d = 3 
provides a correct description of complete wetting in the 
presence of dispersion forces but fails in the regime dom- 
inated by critical fluctuations. This makes it even more 
important to thoroughly analyze a model system which 
takes into account both the bulk and surface critical fluc- 
tuations of semi-infinite systems near T c as well as the 
interfacial fluctuations associated with wetting phenom- 
ena. 

As a representative of the corresponding universality 
class, an Ising ferromagnet is well suited for such an anal- 
ysis. In addition, addressing the issue in two spatial di- 
mensions (d = 2) allows one to gain deep insight because 
the "exact" solvability provides, inter alia, the correct ex- 
ponents. There is a recent experimental interest in such 
genuine two-dimensional systems in the context of pro- 
teins immersed in a fluid two-component lipid membrane 
which is near phase separation, including a critic al p oint 
belonging to the d = 2 Ising universality class jl7l . [T§ | . 
If two such proteins are close to each other the struc- 
tural properties of the membrane between them can be 
described in terms of the strips studied here, giving even 
rise to interesting effective interactions between the pro- 
teins. 

In the present analysis we use the density-matrix renor- 
malization group (DMRG) method 0, Ho, HI, H] to nu- 
merically investigate the behavior of the adsorption in a 
two-dimensional [d = 2) Ising model with short-ranged 
as well as long-ranged surface fields. The DMRG method 
provides essentially exact numerical results for thermo- 
dynamic quantities and correlation functions, including 
the magnetization profiles. This allows us to study sys- 
tematically the adsorption properties along various ther- 
modynamic paths and thus to test the aforementioned 
predictions of the mesoscopic, effective interface Hamil- 
tonian approaches. The DMRG method is based on the 
transfer matrix approach and it is a numerically very ef- 
ficient iterative truncation algorithm for constructing the 
effective transfer matrices for large systems. The method 
was originally developed by White [19|] for the diagonal- 
ization of quantum spin chains. It was then adapted by 
Nishino [20 ] to two-dimensional classical systems. The 
DMRG method allows one to study strips that are in- 
finitely long with widths up to L = 700 lattice constants 
with arbitrary surface and bulk fields. The comparison 
with exact results in the case of a vanishing bulk field 
and in the presence of contact surface fields shows that 
the DMRG method provides a very high accuracy for a 
broad range of temperatures. 

Because the method requires finite values of L the in- 
fluence of the distant surface on the structural properties 
near the surface under study cannot be neglected, in par- 
ticular, near the critical point — even for strips as wide 



as L = 700. On one side, this for the DMRG method 
unavoidable confinement of the system complicates the 
aforementioned crossover behavior due to finite-size ef- 
fects and capillary condensation. On the other side, the 
interplay between complete wetting, critical adsorption, 
and capillary condensation do not appear to have been 
discussed before, although it might be of relevance for 
adsorption-induced colloidal aggregation in binary liquid 
mixtures [HI, [24| . Our approach is naturally suited for 
investigating interesting and relevant aspects of these im- 
portant phenomena. 

The paper is organized as follows. In Sec. [TT] theoreti- 
cal predictions for semi-infinite systems are summarized. 
In Sec. IIIII we describe the microscopic model and the 
method. Numerical results for the phase diagram of the 
present model are described in Sec. IIV1 In Sec. [V] we 
determine the wetting temperature for different ranges 
of the boundary fields. Results for the magnetization 
profiles, the thickness of the wetting layer, and the ad- 
sorption along the isotherms and along the line of the 
pseudo-phase coexistence are reported in Sec. [VI] Sec- 
tion |VlT] summarizes and concludes our results. 

II. DESCRIPTION OF SEMI-INFINITE 

SYSTEMS 

A useful global characterization of the interfacial struc- 
ture near a single wall is provided by the coverage T de- 
fined as the excess number of fluid particles per area ad- 
sorbed on the confining substrate: 

r= / ( P (y)- Pb )dy (3) 
Jo 

where pi, is the bulk number density for a given temper- 
ature T and chemical potential /x; here the fluid number 
density profile p(r) = p(y) is assumed to vary only in the 
direction normal to the wall located at y = 0. 

A. Complete wetting 

If the substrate potential is sufficiently strong there 
is a wetting transition temperature T w such that, if the 
bulk gas phase approaches gas-liquid coexistence po(T) 
along isotherms at temperatures T > T w , T diverges due 
to complete wetting, i.e., a macroscopically thick wet- 
ting film is formed. The equilibrium thickness £q of the 
wetting film can be defined as 

£o = T/( Pl - p g ), (4) 

where pi and p g are the bulk number densities of the 
liquid and gas phase, respectively, at coexistence. In 
magnetic language the fluid is an Ising ferromagnet, the 
gas phase corresponds to the spin down phase, the liq- 
uid phase to the spin up phase, and the difference be- 
tween the substrate potential and its analogue for fluid- 
fluid interactions to a surface field; the undcrsaturation 
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Ap = (p — fio{T))/ (ksTc) is proportional to the bulk 
field H. 

In the complete wetting regime the increase of the ad- 
sorption upon approaching bulk coexistence can be de- 
scribed as Q 

r(A/i-0,T)~ |A/i|-^°, T W <T<T C . (5) 

The exponent fi c ° for this surface quantity depends on 
the form of the fluid-fluid and substrate-fluid forces as 
well as on the spatial dimension d. In d = 3, /5 C ° is 
non-universal; (3f = (i.e., ~ In A/i) for short-ranged 
forces whilst (3 C S ° — 1/p for wall- fluid and fluid- fluid 
pair potentials decaying as r~^ d+p ^ (p = 3 for non- 
retarded dispersion forces). Because the upper critical 
dimension for complete wettin g w ith long-ranged forces 
is d* s = 3 - 4/(p + 1) < 3 Hi| these mean-field 

exponents (3 C S ° = 1/p are not altered by interfacial fluc- 
tuations. For short-ranged forces, i.e., p — > oo one has 
d* s = 3 so that fluctuations matter in <i = 3, but it turns 
that for complete wetting they change only the amplitude 
of the thickness of wetting film [H . 

In d — 2 interfacial fluctuations in the wetting films 
are much stronger and 2 = d < d* for both short-ranged 
and dispersion forces (p = 4 for the latter in d = 2). 
Accordingly the complete wetting exponent takes a uni- 
versal value which turns out to be given by (3f = 1/3 , 
provided p > 3 so that fluctuations dominate [H, [2q . 
From the point of view of an effective interface Hamilto- 
nian this latter universality of (5 C ° is due to the entropic 
effects of the unbinding interface which give rise to an 
effective repulsive interaction for the gas-liquid interface, 
taken to be located on average at y = £, which decays 
- i' K with k = 2(d - l)/(3 - d); n{d = 2) = 2. If 
this entropic repulsion dominates the effective interac- 
tion contribution ~ (,~^~ x \ i.e., if p > k + 1 one finds 
(3 C S ° = 1/3 in d — 2; this defines the so-called weak fluc- 
tuation regime for complete wetting. According to this 
argument, for p < 3, one has (3 C S ° = 1/p. The consider- 
ations leading to the above predictions are valid only if 
the equilibrium wetting film thickness £q is much larger 
than the bulk correlation length £, i.e., £q 3> £• 



B. Critical adsorption 

Near a critical point T c , a confining wall generically 
provides an effective surface field hi acting on the or- 
der parameter (OP) m describing the continuous phase 
transition and leading to the so-called critical adsorp- 
tion P, [3, The ensuing decay of the OP profile 
m(y) = (p(y) — p c )/p c , where p c is the critical density, 
follows the power law 

m(y) ~ jr"/", T = T c , (6) 

where f3 and v are the known critical exponents of the 
bulk OP and of the correlation length £, respectively. Off 



the critical point a crossover to the exponential decay 

m(y) ~ exp(-y/£) (7) 

takes place at the distance y ~ £ from the surface. This 
is described by universal scaling functions P± such that 
m±(y,T,Ap) = m' b a) P±(y/£±(T,Ap),Ap\T\~ A ), where 

to^ ' = m(y — > oo,T < T c ,Ap = 0~) = ?tio|t| /3 < 0, 
t = (T — T c )/T c , and A is the so-called gap exponent; 
± corresponds to r ^ 0. Accordingly the adsorption r± 
diverges as 

r± = m \r\^ ± (r, Ap)K±(Ap\r\~ A ), (8) 

where K±(x) = J™{P±(y±,x)-P±(oo,x)]dy± Near 
the critical point the bulk correlation length acquires the 
scaling form 

e±(r,A M ) = |r|-^ ± (A / i|r|- A ), (9) 

where E± are scaling functions. Since S±(0) = const = 
£q one has 

Z±(t,0)=£\t\- v , (10) 
whereas S±(x — > ±oo) ~ |x| _ly / A and thus at r = 

^(0,Ap) = ^\Ap\-^ A . (11) 

Together with K±(x — ► 0) = const and K±(x — ► ±oo) ~ 
|a;|^' A , this implies that upon approaching (T c , po(T c )) 
along the path Ap = the adsorption diverges according 
to Eq. fl} whereas it diverges according to Eq. © along 
the path r = 0. We note that both Eq. flU and Eq. © 
are consistent with the scaling behavior of T ~ due 
to m ~ |r|' 3 and £ ~ |r| _ " for A/i = and due to 
m ~ lA/il 1 / 5 , (5 = A//3, and £ ~ \Afi\-^ A for r = 0. 
Within MFT j3 = v = 1/2 which results in a logarithmic 
divergence along both paths. 

Equation ^ is expected to hold [l3| also below T c 
for | A/i | 3> |r| A on both the wetting and non- wetting 
side of the coexistence curve, i.e., if the argument of the 
scaling functions E± approaches infinity so that the be- 
havior of £± is governed by Ap (Eq. (fTTj) ). This regime 
is referred to as the critical adsorption regime. Simi- 
larly, above T c Eq. (fTJ) is expected to be valid within 
the regime \Ap\ <C |t| a , i.e., where the behavior of £ is 
governed by r (Eq. (fT0|) ). The same should hold close 
to bulk coexistence below T c . However, it has been ar- 
gued [HI, [ID, Hi] that even near the critical point long- 
ranged forces are important, in the sense that for r fixed 
and A/i — > 0~ the system will always cross over to the 
complete wetting regime (Eq. ((SJ) although upon ap- 
proaching T c (i.e., for smaller r) this crossover will occur 
closer and closer to the coexistence curve. The expecta- 
tion for this scenario to hold is based on the argument 
that long-ranged forces are relevant and thus dominant 
for £ > £. Since £ ~ lA/il"' 3 " and £ - \t\~ v , this 
leads to the condition \ Ap\ < It^/^", i.e., |A/x| < t 3 " 
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for dispersion forces in d = 3. Thus on the gas side 
region between the coexistence curve Ap = and the 
curve A/i = const x |t|"/&> the adsorption should be 
governed by wetting phenomena whereas the divergence 
of the adsorption \t\®~ v according to Eq. JT]) should be 
limited to the range \An\^l A < |r| < \Afi\ P '"/ v . Fur- 
thermore, the effective interface Hamiltonian approach 
predicts that contrary to the critical adsorption regime, 
in the wetting dominated region the divergence should 
depend on the choice of the thermodynamic path taken. 
In particular, for any isotherm within this regime one 
expects T ~ |A/i| _/3 = , but along a path Ap = const\r\ x 
with x > v I {3 C ° one expects T ~ |t|^ -x ^» . This follows 
from Eq. © with Ap ~ \t\ p and £ Q ~ \Ap\-^° . 



C. Crossover phenomena 

These regimes described above are expected to give rise 
to rich crossover phenomena for the adsorption T(Ap, T) 
upon crossing boundary lines between them along various 
thermodynamic paths. Two of them are particularly rel- 
evant for the present work: (i) isotherms T = const < T c 
with A/i — > 0~ (see path (II) in Fig. [T|) and (ii) a path 
parallel to the coexistence curve on the gas side with a 
small undersaturation A/i = const < and with T — > T c 
(path (I) in Fig.©. 

Along an isotherm T = const > T w and below 
the curve A/i ~ — |t| a (see Fig. [T|) T increases as 
| A/i|^ -l/ )/ A upon approaching bulk coexistence from the 
gas side, i.e., for decreasing |A/i|, until one enters the 
crossover region between the lines A/i = — const\r\ A 
and An = ~const\T\ u ^' in which T increases further; 
but therein no specific and well defined power law can 
be expected. Finally, when the path crosses the curve 
A/i = — const\rY^ s it enters the regime governed by 
wetting phenomena so that the adsorption should diverge 
as |A/i| -/3s . For p = p the wall is wet and the adsorp- 
tion is infinite (provided T w <T <T C ). 

The behavior of the adsorption along a path parallel to 
the coexistence curve is expected to be equally rich. Let 
us consider the case that the wetting transition at coex- 
istence is continuous at T = T w . According to the the- 
ory of wetting phenomena Q the adsorption T increases 
smoothly to some finite value (because \Ap\ > 0) upon 
approaching the wetting temperature T w . This increase 
is governed by the scaling laws with respect to An and 
r which are associated with critical wetting. Upon a fur- 
ther increase of the temperature the adsorption should 
slightly decrease because T = ioAp, £q = const due 
to A/x = const, and Ap = pi — p g ~ \t\" until the 
crossover line An = — const\r\ 1 ''^^ to the critical adsorp- 
tion regime is reached. There T should increase ~ |t|' 3 ~ i/ 
until the next crossover line A/i = — const\r\ A is encoun- 
tered beyond which T as a function of t should saturate 
at a certain large value ~ |A/i|^~ l/ )/ A . Finally, after 
passing the right branch of the crossover line An ~ — t a 
(above T c ) the adsorption should decrease ~ r' 3-1 ' upon 



increasing r. 

H 




FIG. 2: Scaling behavior of the capillary condensation 
line H ca (T,L) near bulk criticality (T = T C ,H = 0). For 
r = (T — T c )/T c ^ and L — > oo, H ca {r, L) approaches bulk 
coexistence H — as |r| A_I/ /L and the capillary critical point 
(T C}L ,H C , L ) approaches T c and H = as L~ 1/v and L~ A/l/ , 
respectively. The latter is valid for strong surface field and 
takes the form ^ _ ( A_A i)/ I/ i where Ai is the surface gap ex- 
ponent, for weak hi [4(|. Here T Ci l < T c . Note the upward 
bent of the capillary condensation line; this feature is due to 
A — v < 1 (see, c.f., the text following Eq. (|25|) h However, this 
bent is weak: A - v = 7/8 = 0.875 in d = 2 and A - v ~ 0.936 
in d = 3. 



III. MICROSCOPIC MODEL 

In this section we introduce the microscopic model 
within which we investigate quantitatively the crossover 
regimes described above. Specifically, we consider an 
Ising ferromagnet in a slit geometry subject to identi- 
cal boundary fields. Contingent on the type of numerical 
approach we shall use, our results refer to d = 2 strips 
defined on a square lattice of size M x L, M — > oo. The 
lattice consists of L parallel rows at spacing a, so that 
the width of the strip is La; in the following we set a = 1. 
Successive rows are labeled by an index j. At each site, 
labeled (k,j), there is an Ising spin variable taking the 
value o~k.j = ±1. The boundary surfaces are located in 
the rows j = 1 and j — L and periodic boundary condi- 
tions (PBCs) are assumed in the lateral x direction. Our 
model Hamiltonian for the strip with PBCs and AI — > oo 
is given by 



+ E^£^+ F 5>H > ( 12 ) 

3 = 1 k k,j J 

where the first sum is over all nearest-neighbor pairs and 
the external potential is measured in units of J > 0. 
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V^j?* = + V£ +1 _j is the total boundary field experi- 
enced by a spin in the row j; it is the sum of the two in- 
dependent wall contributions. The single-boundary field 
V? is assumed to have the form 



with p > and h\ > 0. H is a bulk magnetic field. 
According to Eq. (Ti"2|) /ii and H are dimensionless. 

This model can be viewed as being obtained from a 
2d lattice gas model mimicking a two-dimensional one- 
component fluid with a short-ranged interaction poten- 
tial between the fluid particles and either short-ranged 
or long-ranged substrate potentials. The equivalence be- 
tween the lattice gas and the Ising model implies the 
following relationships (see, e.g., Ref. [2!|): the bulk 
magnetic field H in the former is proportional to the 
deviation of the chemical potential from the bulk phase 
boundary po(T) in the latter, i.e., H ~ A/i. The lattice- 
gas analogue of the number density in the fluid is re- 
lated to the magnetization by p = (m + l)/2, so that 
Ap = pi — p g ~ 2m&, where is the spontaneous mag- 
netization. Finally 4J corresponds to the strength of 
the attractive pair potential between the fluid particles, 
taken to be short-ranged so that in the lattice gas model 
it can be modelled by nearest neighbor interactions. V ext 
is a combination of the substrate potential and the liquid- 
liquid interaction. These relationships can be extended 
to binary liquid mixtures (30j . 

We recall that although there is no longer any true 
phase transition for finite L, in two-dimensional Ising 
strips with large L there is still a line of extremely weakly 
rounded first-order transitions ending at a pscudocritical 
point the location of which in the plane (H, T) spanned 
by the bulk field and the temperature T depends on the 
character of the surface fields |3ll . [HJ . For surfaces which 
prefer the same bulk phase, the phenomenon equivalent 
to capillary condensation takes place. The pseudo-phase 
coexistence between phases of spin up and spin down oc- 
curs along a line H ca (T, L), which is given approximately 
by the analogue of the Kelvin equation [33j : 

H ca (T,L) « -a(T)/ (L\m<?\T)\) , (14) 

where a(T) is the interfacial tension (divided by J) be- 
tween the coexisting bulk phases and m£ (T) < is 
the spontaneous bulk magnetization for H = 0" . For 
the d = 2 Ising model, the surface tension is given ex- 
actly f3Ja{T) = 2{K — K*), where K = J/k B T and K* 
satisfies sinh 2K sinh 2K* = 1. The occurrence of thick 
wetting films of + spins at the two surfaces, for a ther- 
modynamic bulk state corresponding to — spins, gives 
rise to nontrivial corrections to Eq. (frjj which shift the 
condensation line to larger values of \H\ HU. The 
pseudo-coexistence line ends at a pseudo-capillary crit- 
ical point (-ff c ,L, ?c,l) where T c l(/ii,p) lies below the 
temperature T c of the bulk critical point. Its position as 



well as that of H c ,l depends on L and the strength and 
the range of the surface fields. For large L and strong h\ 
the shifts of the critical temperature and of the bulk field 
H c ,h{hi,p) are given by [H,!!! (see, c.f., Fig. [2) 

T c .l -T c ~ -L- X l\ H C . L ~ -L- A /" (15) 

with v = 1 and A = 15/8 for the two-dimensional Ising 
model. 

IV. PHASE DIAGRAM 

In Fig. [3] we show the phase diagram for the present 
model calculated by using the DMRG method for a strip 
of width 340 with h\ = 0.8 and three choices of the pa- 
rameter p describing the decay of the boundary field: 
p = 50, 3, and 2. The case p = 50 is expected to resem- 
ble the behavior corresponding to short-ranged surface 
forces. In this figure we display the various crossover 
lines discussed in Sects. HI and UT1 and the thermodynamic 
paths along which we have calculated the adsorption F. 
The pseudo-phase coexistence line, the crossover lines, 
and the adsorption have been determined for 2d Ising 
strips within the DMRG method. 

The thick solid line in Fig. [3] indicates the bulk phase 
coexistence line (if = 0, T < T c ) terminating at the bulk 
critical point (H = 0,T = T c = [ln(l + y/2)]- 1 J/k B ~ 
2.269J//cb)- The symbols (open circles for p = 2, open 
squares for p = 3, and stars for p = 50) show the pseudo- 
phase coexistence. It turns out that the pseudo-phase co- 
existence line for long-ranged boundary fields is located 
slightly further away from the bulk coexistence line, espe- 
cially at lower temperatures, than the pseudo-phase coex- 
istence line for the short-ranged boundary fields (p = 50). 
These pseudo-lines have been identified as those posi- 
tions (H, T) in the phase diagram where the total mag- 
netization of the strip vanishes, i.e., X)j=i 171 j = with 
rrij = Kj). 

At fixed T (fixed H) at those positions the free en- 
ergy of the strip exhibits a maximum as a function of H 
(T) which is the rounded remnant of the nonanalyticity 
(kink) of the free energy at the first-order phase tran- 
sition in the bulk [L = oo) system. As can be seen in 
Fig. H for T <T C (T > T c ) the line defined by the zeros 
of the total magnetization moves to less (more) nega- 
tive values of H upon increasing T. Determining the 
pscudocritical temperature T Cj l, where the pseudo-phase 
coexistence line ends, is difficult because in this quasi- 
onc-dimcnsional system the critical point is not a sharp 
concept and one has to examine various criteria in order 
to estimate T Cj l- However, these different criteria, such 
as the maximum of the specific heat, or the erosion of the 
jumps in the adsorption and in the solvation force upon 
crossing the line H ca (T, L;p), provide somewhat different 
estimates for T C] £. Determining T c .l for different values 
of p and L is beyond the scope of the present analysis. 
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FIG. 3: Phase diagram for a d = 2 Ising strip subject to iden- 
tical boundary fields Vf xt (see Eq. |T3|) obtained by using 
DMRG for a strip width L — 340 and the amplitude hi = 0.8 
for the boundary fields. The thick solid line indicates the 
bulk coexistence line. The lines interpolating the symbols 
represent the pseudo-phase coexistence lines H ca (T, L;p) for 
three different values of the exponent p governing the alge- 
braic decay of the boundary fields: open circles correspond 
to p = 2, open squares correspond to p = 3, and stars cor- 
respond to p = 50. Various crossover lines and 
with di = 1 = d,2 = 1 and v = 1 and A = 15/8 in d = 2, 
as discussed in the main text, are shown as well as the ther- 
modynamic paths (0) — (4) along which the adsorption has 
been calculated. (1) — (4) are various isotherms and (0) runs 
parallel to the pseudo-coexistence line. Thus path (0) resem- 
bles (in a certain sense, see later) a path of type (I) in Fig. [1] 
whereas the paths (1) — (4) do correspond to a path of type 
(II) in Fig. □ 



For the width L = 340 and hi = 0.8 we approximately 
obtain T Cj l w 2.2; for more details concerning the short- 
ranged case see Ref. [36[. 

We have investigated how the location of the pseudo- 
phase coexistence line iJ CQ (T, L;p) of capillary conden- 
sation changes as a function of p for fixed width L, 
fixed amplitude hi of the surface field, and fixed tem- 
perature T. In Fig. 0] we plot the logarithm of the 
difference H ca (T, L;p) — H ca (T, L;oo) as a function of 
p calculated for hi — 0.8 and (i) for several values of 
the temperature at fixed L = 300 (Fig. Ufa)) as well 
as (ii) for several values of L at fixed T = lSJ/ks 
(Fig. HJb)). (Note that p = oo corresponds to a pure sur- 
face contact field at j = 1, L (see Eq. (fl~3|) ).) From these 
plots we can clearly distinguish two different regimes 
for the behavior of the pseudo-phase coexistence line, 
which are separated by a crossover region occurring for 
1 ~ P ~ 3. In both of these two regimes the shift 
H ca (T, L;p) — H ca (T, L; oo) varies exponentially but with 
different decay constants. Moreover, for p < 1 the shift 
relative to the short-ranged pseudo-phase coexistence line 
does not depend on the temperature; the dependence 
on the strip width L becomes, for the range of L con- 
sidered here, very weak and finally negligible for p < 



with A ~ 5.2 and Bq ~ 1.5. For p > 3 one has 
H ca (T,L;p) - H ca (T,L; oo) ~ B(L, T) exp(-Cp) with 
C ~ 0.723. We note that, as expected, for p fixed wider 
strips and higher temperatures give rise to smaller shifts 
of the pseudo-phase coexistence line. 
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FIG. 4: The relative shift of the pseudo-coexistence line for 
Ising strips as function of the range of the boundary fields 
characterized by the decay exponent p (Eq. (|13[0 . The cal- 
culations are carried out by using the DMRG method for 
hi = 0.8 and (a) at fixed L = 300 and for four temperatures 
T* = k B T/J; (b) at fixed temperature (k B T/J = T* = 1.8) 
and for four thicknesses L of the strip. 



V. DETERMINATION OF THE WETTING 
TEMPERATURE FOR DIFFERENT RANGES OF 
THE BOUNDARY FIELDS 

It is important to assure that the range of tempera- 
tures, for which we perform our calculations of the ad- 
sorption, lies above the wetting temperature T w (jp,hi). 
For hi = 0.8 and a short range of the surface fields (p = 
oo) the wetting temperature in the semi-infinite d = 2 
Ising system equals T w (p = oo, hi = 0.8) ~ 1.41J/fcg 
and the wetting transition is second order. Adding the 
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long-ranged tail to the boundary fields shifts the wet- 
ting temperature. Wetting transitions in d = 2 in the 
presence of long-ranged substrate potentials have been 
discussed by Kroll and Lipowsky [37j within the contin- 
uum one-dimensional solid-on-solid model, with the con- 
clusion that there is no wetting transition for fields which 
decay more slowly than 1/j 3 . In this case their predic- 
tion is that the interface remains pinned to the wall at all 
finite temperatures. This conclusion has been supported 
by providing upper and lower bounds for the ground- 
state energy E of the corresponding Schrodinger equa- 
tion. Within this approach the existence of the bound- 
state solution of the Schrodinger equation (E < 0) 
corresponds to a localized interface. When the bound 
state ceases to exists, Eq — > signals the critical wet- 
ting transition. The construction of upper and lower 
bounds (based on using an exponentially decaying trial 
wave function) shows that the ground-state energy has 
a non-zero value (i.e., there is no transition) for a finite 
potential strength. The case p = 2 has been studied by 
Privman and Svrakic in Ref. [H[ also within the SOS 
model; they argued that for attractive effective poten- 
tials decaying like 1/r the wetting transition is no longer 
sharp (becomes rounded). For the d — 2 Ising fcrromag- 
net the wetting phenomena for the case of the marginal 
range p = 3 have been studied within the DMRG ap- 
proach Hi] and by MC simulations [10, [H|. These stud- 
ies have shown that the presence of long-ranged tails in 
the boundary fields decreases the wetting temperature 
with respect to the case of the short-ranged surface fields 
and that T w (p = 3, hi = 0.8) ~ 0.75J/&B- Apart from 
Refs. [37l. |38|. the case p < 3 has not been studied so far. 
In order to check the predictions by Kroll and 



Lipowsky, we investigate whether in the d = 2 Ising fer- 
romagnet with p = 2 a wetting transition exists. In the 
finite systems we are studying, below the bulk critical 
temperature the wetting transition as a function of hi 
for any range p of the boundary fields can be inferred 
from the so-called (weakly rounded in d = 2) interface 
localization-delocalization (ILD) transition [H,[4^ which 
occurs in strips with antisymmetric boundary fields, i.e., 
for Vf xt = Vf — V£ + i_j. This transition is the precur- 
sor of a wetting phase transition that occurs in the limit 
of infinite film thickness (L — > 00) at the critical curve 
T w {p,h\). For T < T w (p,hi)(T > T w {p,h\)) such an 
interface is bound to (unbound from) the walls [HI, |42j . 

First, we note that T w {h\ — > 0,p) — > T c for any p (see 
Eq. H2])). On the other hand hi w (T = 0,p), which is 
the solution of the implicit equation T w (h\,p) = and 
denotes the critical surface field strength beyond which 
the system is wet even at T = 0, shifts towards lower 
values upon decreasing p. There is no reason to expect 
a non-monotonic behavior of T w {h\,p). Therefore, the 
gross features of the shape of the wetting transition line 
T w (h-\_,p) for an arbitrary p > 1 can be inferred from 
localizing the position of the wetting transition h\ w (T = 
0,p) in the ground state. The ground-state energy of 
the system (in units of J) can be found directly from 
the Hamiltonian with a vanishing bulk field. Because 
the system is translationally invariant along a strip, it is 
sufficient to consider only the configurations of a single 
column. 

In a partial wetting regime with antisymmetric surface 
fields there are only two coexisting states (all spins up or 
all spins down) with the energy (per the number M of 
columns) 
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-J(L-l)ThiJ [ 1- 
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LP 2p (L - 1)P 



(L - 1)p 2p 



1 

L~p 



= -J(L-l). 



(16) 



In a wet regime at least one interface between spin up and 
spin down configurations has to be present. For even L 
the lowest-energy configuration is that for the state with 
the interface located in the middle of the strip. For p — > 
00 (short-ranged case) L — 1 degenerate states emerge 



with a single interface positioned at any of the rows but 
the ones closest to the surfaces. In general, the energy 
of a single column with an interface in the middle of 
a strip (spin up for j = 1, . . . , L/2 and spin down for 
J /• 2 • ! /.:■ is 



E wet = -J(L - 3) - hi J I 1 - -L + 1 



LP 2P {L-1)p "' {L/2)p {L/2 + \)p 



{L/2 + l)P (L/2)p "' {L~1)p 2P LP 

I 

In the ground state and at the ILD transition the en- crgies E pw and E wet are equal. For a particular L this 
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determines the magnitude (£): 



h[^(L) =1/ (! + _ + _ + ...+ 



1 

3^" 



(L/2)p (L/2 + l)P 



(L - 1)p LP 



oo 1 oo oo 1 

E-- 2 E -+ E - 



n=£/2+l 



(18) 
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In order to find the critical wetting field h\ w (L — ► oo) we 
take the limit L — > oo. According to the second line in 
Eq. (JTHK this leads to 



OO 

1 /E^ = vc(p) 



(19) 



n=0 



where £(p) is the Riemann zeta function. Its values are 
known analytically only for certain even values: £(p = 
2) = tt 2 /6, ((p = 4) = tt 4 /90 which gives h lw (p = 2) « 
0.6079 and /n w (p = 4) w 0.9239. Other values can be 
found in tables of special functions, e.g., Q{p = 3) w 
1.2021 giving h\ w (p = 3) r* 0.8319. In the short-ranged 
limit (p — > oo) the Riemann zeta function approaches 1, 
confirming Abraham's solution at T = [43j . Even more 
interesting is the opposite limit £(p — > 1) — > oo which 
results in hi w (p — > 1) — ► 0. This means that at T = 
the wetting transition does not exist for p < 1. 
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FIG. 5: Phase diagram at bulk coexistence H = for con- 
tinuous wetting transitions in the d = 2 Ising ferromagnet 
for (Eq. I|13| l) p = 2, p = 3, and p = 50 (quasi-short-ranged 
boundary field) obtained within the DMRG approach. For 
p > 1 the values in the limit T — » are known exactly 
for semi-infinite systems. The dashed line is the analytically 
known exact result for p = oo and for semi-infinite systems 
0. T* = k B T w /J. 

In order to determine the location of a quasi-ILD tran- 



sition between T = and T = T c for p < 3 various 
criteria can be applied [3^]. Here we have adopted the 
approach involving the magnetic susceptibility \- The 
singularity (or a maximum) of the magnetic susceptibil- 
ity x is one of the most useful criteria for the localization 
of a phase transition (or of a pseudo-phase transition for 
a finite system). The magnetic susceptibility can be cal- 
culated as the second derivative of the free energy / with 
respect to the bulk magnetic field H . This method is very 
convenient for the DMRG approach because the latter 
provides the free energy with a very high accuracy. Nev- 
ertheless the present case is somewhat special, because 
we want to determine T w (p, hi) at H = (see Fig. [5]), 
where in the partial wetting regime, i.e., for T < T w , 
there is a first-order bulk transition. In the thermody- 
namic limit there is coexistence of phases with opposite 
magnetizations. Thus there is a discontinuity of the first 
derivative of the free energy / (a jump of the magnetiza- 
tion m = —df/dH) upon changing the sign of the bulk 
magnetic field. Accordingly, in order to calculate x there, 
one has to calculate the derivatives for small nonzero bulk 
fields and then to consider the limit H — ► 0. In the com- 
plete wetting regime, i.e., T > T w , or equivalently above 
the ILD transition, the finite system exhibits only a sin- 
gle phase with an interface meandering freely between 
the walls so that there is no discontinuity of the free en- 
ergy derivatives upon crossing H = 0. For numerical 
calculations such as the ones within the DMRG method 
the necessity of performing an extra limiting procedure 
(H — > in this case) is cumbersome. Therefore, instead 
of x, we have focused on another quantity xoi which also 
corresponds to the second derivative of the free energy at 
fixed T and hi, but is calculated numerically in a sym- 
metrical way with respect to H = by taking the free 
energy values at five equidistance points: — 2AH, —AH, 
0, AH or 2AH; we typically used AH = 10 5 . Be- 
cause our calculations are always carried out for finite L, 
there is no discontinuity of the magnetization in the par- 
tial wetting regime. These discontinuities are replaced 
by functions which are rounded but steeply varying at 
H = 0. In order to determine the ILD transition we 
have scanned the phase diagram at fixed hi. The higher 
the temperature, the less steeply the magnetizations vary 
and the values of their derivative xo are smaller. Above 
the wetting temperature, where there is no discontinuity, 
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Xo{H = 0, T) saturates for increasing T; here xo is equiv- 
alent to x- Therefore, at fixed L, the ILD transition can 
be identified by the maximal slope of xa or the minimum 
of its derivative with respect to temperature. Although 
all derivatives have been performed numerically, the high 
accuracy of the DMRG method guarantees very precise 
results. 
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FIG. 6: Magnetization profiles, relative to their bulk values 
rrib(T,H) < 0, near one wall (j > 1) along the isotherm 
T* = 1.8 (see path (1) in Fig.E}, calculated within the DMRG 
method for d = 2 Ising strips of width L — 500, for (a) short- 
ranged (p = 50) and (b) for long-ranged (p — 2) boundary 
fields of strength hi — 0.8. Along this isotherm the pseudo- 
capillary condensation transition occurs at H ca — —0.00172 
for p = 50 and at H ca ~ -0.00195(4) for p = 2. Thus both in 
(a) and (b) the full line for H = -0.0017 and H = -0.0019, 
respectively, corresponds to the capillary filled state with 
m(j)-m b (T,H)~ 2\m b (T,H)\. 



Finally we have extrapolated T w {p 1 hi, L) to the limit 
L — > oo in order to obtain the wetting temperatures 
T w (p,hi). The obtained wetting phase diagram in the 
(T, hi) plane is shown in Fig. [5] for p = 2, 3, and p = 50. 
The dashed line is the exact result by Abraham for p = oo 
[IK . The close agreement of our data for p — 50 with this 
exact result generates confidence in our numerical proce- 
dures. Our results show that contrary to the claim by 



Kroll and Lipowsky [37[ also for p = 2 there are wetting 
transitions and that the wetting temperature of these 
continuous transitions can be well localized as a function 
of hi. In view of the calculations of the adsorption, for 
our choice hi = 0.8 the surface is wetted for all temper- 
atures. As already mentioned above, the claim by Kroll 
and Lipowsky is based on using a continuum planar solid- 
on-solid (SOS) approximation. SOS models ignore bulk- 
like fluctuations such as bubbles and the formation of 
" overhangs" of the line separating oppositely magnetized 
regions, while focusing on describing the long- wavelength 
behavior of the interface, expected to play the crucial role 
in wetting phenomena. In the present case of the d = 2 
Ising model the bulk fluctuations arc particularly strong 
giving rise to a very diffuse interface. For the case of 
p = 2 they might provide the mechanism for the unbind- 
ing of the interface from the wall, even though the long- 
wavelength interfacial fluctuations are strongly supressed 
by this very long-ranged boundary field. 



VI. NUMERICAL RESULTS AND DISCUSSION 

Using the DMRG method we have performed calcula- 
tions of the magnetization profiles m(j) from which we 
have obtained the thicknesses £o of the wetting layer and 
the adsorption T. The analysis of the shapes of the pro- 
files m(j) and of the scaling properties of £o gives a better 
understanding of the behavior of T along different ther- 
modynamic paths. We recall, that the predictions for 
the scaling behavior of the adsorption follows from the 
scaling behavior of £q. 

In order to infer possible power laws governing the be- 
havior of the thickness of the wetting layer and of the 
adsorption along different paths we have calculated local 
exponents of the quantities of interest as a function of H 
or t. They are defined as 



lnQ(» + l) -lnQ(i) 



lnxi+i - lnxi 



(20) 



which is the discrete derivative of data Q as a function 
of x in a log- log plot; here Q = £q or T and x = H 
or t. Such a quantity probes the local slope at a given 
value Xi of x at a point i along the path considered. It 
provides a better estimate of the leading exponent than 
a log-log plot itself. We have chosen the convention that 
if Q decays algebraically as a function of x then Q ~ x~ z 
with z > 0. The high quality of the DMRG data allows 
us to reliably carry out this numerical differentiation. 



A. Isotherms 

1. Magnetization profiles 

In Figs. [5] and [7] we show a selection of magnetization 
profiles calculated for a strip width L = 500 with p = 50 
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and 2 along the two isotherms T* = 1.8 and T* = 2.25 
indicated in Fig. [3] as path (1) and (3), respectively. 
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FIG. 7: Magnetization profiles, relative to their bulk values 
rrib(T,H) < 0, near one wall along the isotherm T* = 2.25 
(see path (3) in Fig. [3} , calculated within the DMRG method 
for d = 2 Ising strips of width L = 500, (a) for short-ranged 
(p — 50) and (b) for long-ranged (p = 2) boundary fields of 
strength hi — 0.8. Along this isotherm the pseudo-capillary 
condensation transition occurs at H ca — —0.000191 for p — 50 
and at H ca ~ -0.00026(4) for p = 2. Thus both in (a) and 
(b) the full line for H — —0.0001 corresponds to the capillary 
filled state with m(j) - m b (T,H) ~ 2\m b (T,H)\. 

From these plots one can see that within the acces- 
sible range of values H < H ca {T 7 L;p) the shapes of 
the profiles are not slab-like, even though the system 
size is very big so that \H ca (T,L)\ is expected to be 
sufficiently small. They are characterized by a pro- 
nounced spatial variation at the emerging interface be- 
tween the spin up and the spin down phase, where the 
profile varies quasi linearly with the distance from the 
wall, and by extended tails. For the short-ranged case 
(p = 50) the profiles along the low-temperature isotherm 
(T* = 1.8) exhibit a rapid decay to their bulk values 
m b = m b {T,H) (< for H < 0) with the bulk correla- 
tion length $(T,H) = |r|-"5:(H|r| A ) (see Eqs. ©-G2J) 
as decay length which increases upon approaching the 
pseudo-phase coexistence. Only very close to the pseudo- 
phase coexistence line H ca (T, L;p) the emerging shape of 



the profile to a certain extent resembles that characteris- 
tic of the free interface between the spin up and the spin 
down bulk phases with a broad interfacial region and 
a rather narrow region where the magnetization stays 
somewhat close to the value —m b > characteristic for 
the wetting phase corresponding to H = + ; this latter 
region thickens as H — > H ca (T, L;p). 

For the long-ranged forces (p — 2) the profiles along 
the low-tempcraturc isotherm (T* — 1.8) exhibit an 
interface-like shape already away from the pseudo-phase 
coexistence line H ca (T, L;p). The interfacial part of these 
profiles is much more pronounced than for the short- 
ranged case (p = 50) and their decay is much slower. The 
exponential decay of the profiles to the value rrib(T,H) 
ultimately crosses over to the power law decay ~ j~ 2 . In 
general, algebraically decaying inter-particle and surface 
fields are known to generate algebraically decaying order 
parameter profiles [44| . 

For both p = 50 and p = 2 the profiles corresponding to 
H — > along the high-temperature isotherm (T* = 2.25) 
decay much slower as compared to the low-temperature 
isotherm (T* = 1.8) and their shapes become quantita- 
tively different (see Fig. [7]). The magnetization in the 
first few layers near the walls decreases more rapidly so 
that there the profiles acquire a positive curvature. Fur- 
ther away from the wall there is an inflection point and 
the narrow pleateau occurring near the boundaries for 
the low-temperature isotherm (see Fig.|6jb)) disappears. 
The interfacial region becomes very broad. 



2. Thickness of the wetting layer 

In order to infer equilibrium thicknesses £q of the wet- 
ting layers from the magnetization profiles we have cho- 
sen the criterion of the change of sign of the curvature of 
the profile, i.e., we have assigned a thickness to a wetting 
layer which corresponds to the distance j = £q at which 
the profile exhibits its inflection point. Adopting other 
criteria for determining the wetting film thickness, e.g., 
the thickness which corresponds to the distance j = jo 
at which the magnetization vanishes, m(jo) = 0, leads, 
as it must be, to the same conclusions about physical 
observables. 

We have studied the behavior of £q along several 
isotherms. Representative data are shown and discussed 
below. They correspond to the isotherms T* = 1.8,2.25 
(see the paths (1) and (3), respectively, in Fig. [3]) and 
T = T c (see path (4) in Fig. [3]) and have been obtained for 
several strip widths and ranges p = 50, 5, 4, 3, 2, 1.5 of the 
boundary fields. (The data for the isotherms T* = 1.6, 
1.9, 2, and 2.2 are not shown.) We have chosen hi = 0.8 
for which T*,(p = oo) ~ 1.41. All considered isotherms 
lie above the corresponding wetting temperatures T w (p). 

Along the isotherm T* = 1.6 the wetting layers are so 
thin that it is difficult to assign a suitable thickness £o, 
especially for p = 50. Along the isotherm at the higher 
temperature T* = 1.8 (path (1) in Fig. [3]) this assignment 
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is much clearer and the results for £o(H) and its local 
exponents are shown in Fig. [5J They were obtained for 
strips of width L = 500 with h\ = 0.8 and for various 
ranges p of the boundary fields. 
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FIG. 8: (a) Equilibrium wetting layer thickness £q in units 
of the lattice constant a as a function of the bulk magnetic 
field H calculated along the isotherm T* = 1.8 for L = 500, 
h\ — 0.8 and for the ranges p = 1.5,2,3,4, and 50 of the 
boundary fields. A rapid increase of to upon approaching the 
pseudo-phase coexistence line H ca (T,L;p) supports the exis- 
tence of the wetting transition for p = 2 and 1.5. Fhe values 
of H ca (T, L;p) are given by the points at which the various 
lines end on the right. The inset shows a log- log plot of these 
data. The thin solid lines indicate slopes (from top to bot- 
tom): 2/3, 1/2, and 1/3. (b) The local exponents z for £o(H) 
(see Eq. (|20p ). For H — * 0, without the occurrence of capil- 
lary condensation, i.e., for L — > oo, the expected exponents 
are f5 c a ° = 1/3 (*) for p = 3, 4, and 50, 1/2 (■) for p = 2, 
and 2/3 (•) for p = 1.5. The dotted lines indicate the ex- 
pected extrapolations H — * for L — > oo. The log-log plot 
supports these conclusions. We note that the local exponents 
do not always follow from naive visual impressions. For ex- 
ample, z(p — 50) > z(p = 2,3,4), although in (a) the curves 
of £o(p = 2, 3, 4) seem to increase somewhat stronger than the 
one for £o (p — 50) . 

Path (1) lies entirely within the region between the 
bulk coexistence curve and the crossover boundary 



= -di|r| 3l/ with v = 1 for the d = 2 Ising 
model, i.e., within the complete wetting regime (see 
Fig. Indeed, for all values of p we observe the in- 
crease of £q upon approaching the pseudo-phase coexis- 
tence line H ca (T, L;p). The lines in Figs. [8lfl3l end on 
the right when they reach the pseudo-phase coexistence 
line H ca (T, L;p) (see Fig. [3]). Apart from this size depen- 
dence of the iJ-window with access for complete wetting, 
for L = 500 and 300 (data for L = 300 are not shown) 
the finite-size effects are negligible. As expected from 
the shape of the profiles (see Figs. [5] and [7]) for the same 
value of H the wetting layer is much thicker for small 
values of the exponent p than for short-ranged surface 
fields (p = 50). The rate of increase of £o, as character- 
ized by the behavior of the local exponents z, depends 
on the range of the boundary field p and varies along 
the isotherm. For p = 1.5 the rate of increase is al- 
ways largest, i.e., z remains the largest along the entire 
isotherm (see Fig. Ob)). For p = 50,4, and 3 the local 
exponent z decreases as \H\ — > 0, i.e., the initial increase 
of the wetting layer is larger than the one observed very 
close to H = for which the plateau in the profile near 
the wall corresponding to the wetting phase at H = + 
starts to form. For p = 2 the very long-ranged nature of 
the boundary field enforces the occurrence of a narrow 
plateau in the profile near the wall even far away from 
H = 0. Upon decreasing \H\ the width of this plateau 
increases and the emerging interfacial region broadens 
giving rise to the increase of £q which is initially smaller 
than for p = 50. Only very close to H = this trend re- 
verses and eventually the local exponents z corresponding 
to p = 2 become larger. For p = 3 and 4 the rate of in- 
crease of £o remains smaller than for p = 2 and 50 along 
the whole isotherm (for \H\ < 0.008). According to the 
predictions summarized in Sec. |TT] for semi-infinite sys- 
tems (L — > oo), upon approaching the bulk coexistence 
line H = the thickness of the wetting film is expected 
to diverge for p > 3 with the exponent (3f = 1/3 and 
for p < 3 with the exponent (3f = 1/p. As can be in- 
ferred from Fig. [5] the local exponents z calculated for 
p = 50, 4, 3, and 2 tend towards their predicted values 
for H — > but cannot reach them due to capillary con- 
densation. For p = 1.5 the local exponents have to be- 
have non-monotonically in order to reach the expected 
value 1/p = 2/3. The behavior of £q along the isotherm 
T* = 1.9 is very similar (data not shown). 

Path (2) corresponds to T* = 2.0 and runs in be- 
tween the crossover lines given by H£/ = — di\r\ 31 ' and 

He? = -d 2 |r| A with v = 1 and A = 15/8 in d = 2; 
it hits the pseudo-phase coexistence line of capillary 
condensation before it reaches the crossover boundary 
Har = — cZi |t| 3 ■ Along that path the variation of £o 
is very similar to the one observed along the isotherm 
T* = 1.8 with the local exponents approaching the ex- 
pected values of the complete wetting regime (data not 
shown). Again the asymptotic behavior is preempted by 
capillary condensation. 
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FIG. 9: (a) Wetting layer thickness £o in units of the lattice 
constant a as a function of the bulk magnetic field H calcu- 
lated along the isotherm T* = 2.25 (path (3) in Fig. [3| for 
L = 500, hi = 0.8, and for the ranges p = 1.5, 2, 3, 4, and 50 of 
the boundary fields. On the right the curves end at the corre- 
sponding capillary condensation point H ca (T, L;p). The inset 
shows a log-log plot of these data. The thin solid line indi- 
cates the slope 8/15. (b) The local exponents z for £o(H) (see 
Eq. J20J). For p > 2 and H -> they tend towards 0.48, 
0.52 whereas the expected value 8/15 ( k) is somewhat higher. 
According to Fig. [3] along this isotherm capillary conden- 
sation prevents one to reach the ultimate complete wetting 
regime. For p = 1.5 the behaviors of £o(H) and z(H —> 0) are 
distinctly different from those for p > 2. 

Results along the near-critical (T* = 2.25) and critical 
(T* = T*) isotherms (paths (3) and (4) in Fig. |3J) are 
displayed in Fig. [5] and Fig.[Tni respectively, for L = 500. 
(For L = 300 we observe a similar behavior; these data 
are not shown.) Both paths lie in the critical adsorption 
regime where one expects lo{H — ■> 0) ~ £ ~ \AH\~ U / A 
with u/A = 8/15 in d = 2. This universal behavior is 
expected to hold if the long-ranged boundary field decays 
sufficiently rapidly, i.e., if p > (d + 2 — rj)/2. In this case 
the long-ranged part of the boundary field is irrelevant 
in the RG sense with respect to a pure contact surface 
field Here r\ is the critical exponent governing the 
algebraic decay of the two-point correlation function in 
the bulk and at T c . In the present case of the d = 2 Ising 
model rj(d = 2) = 1/4 so that for p > 15/8 we expect to 



observe the power law £q ~ |AiJ|~ 8 / 15 . 

For p = 1.5 one has p < (d + 2 — rj)/2. To our knowl- 
edge so far this case has not been studied in the literature 
in the context of critical adsorption. According to gen- 
eral scaling arguments the magnetization profile for 
the semi-infinite system under rescaling of distances by 
a factor b behave as 

m±{y,T,H,h 1 ,p) = 

b- f) /"m±{yb- 1 ,Tb l l v ,Hb A ' v Mb" 8 ), 

(21) 

where ui s = (d + 2 — rf)/2 — p. Because in the present 
case u> 8 > the range of such wall-fluid interactions is 
relevant in the RG sense at and near T c , i.e., the profiles 
are affected by the decay of the boundary field even if y 
is sufficiently large; in the critical region the profiles are 
no longer equivalent to those generated by a pure contact 
surface field (p = oo). Choosing b ~ y one obtains the 
scaling form of the magnetization 

m±(l/,r,ff,/n,p) = y- /3/v M±(y/Z T ,y/Z H ,y"'h 1 ), 

(22) 

where £ r = £±(t, 0) is given by Eq. (JTUJ) whereas £h = 
£(0, Afi) is given by Eq. (fTTj) . For oj s < and large 
enough y the scaling function A4± can be expanded in 
terms of their third argument demonstrating that in this 
case the long range of the wall-fluid interaction gives rise 
to contributions to the scaling functions of the magneti- 
zation profile which are subleading to the dominant uni- 
versal y~^l v behavior valid for y < £ and p = oo. Such 
an expansion cannot be performed for large y if lo s > 
and in order to obtain a prediction in this case for the 
leading behavior of the magnetization profile and hence 
the thickness of the wetting layer £$ and the adsorption, 
one has to calculate the full generalized scaling functions 
of the magnetization profiles. For p < 15/8 such an anal- 
ysis represents a necessary future research goal which, 
however, is beyond the scope of the present study. Never- 
theless, we find it instructive and stimulating to show our 
numerical results for p = 1.5 because they are strikingly 
different from those obtained for values of p satisfying 
p> {d + 2-rj)/2. 

As along the noncritical isotherms T* < 2.25 we find 
that also along the isotherms T* = 2.25 and T = T* the 
increase of the wetting film thickness upon approaching 
H = is much more stronger for p = 1.5 than for the 
other values of p. For \H\ < 0.01 the local exponents 
attain z(p = 1.5) ~ 0.64, which indicates an algebraic 
increase of £o in this range of H with £q f=a |iJ|~°- 64 . 
This effective exponent is close to the value 1/p = 2/3 = 
0.66(6) predicted for the complete wetting behavior for 
p = 1.5. Closer to H = the functional form of the 
increase changes and z increases strongly. 

For p = 50,4, and 3 the local exponents decrease 
upon approaching the pseudo-phase coexistence line 
H ca (T, L;p). This behavior is similar to the one occur- 
ring along the isotherm T* = 1.8; however, apart from 
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the case p = 1.5 the values of z in Fig. [Hb) arc larger 
than the values in Fig.[5{b), which indicates the stronger 
increase of the wetting film thickness. For p — 2 the 
variation of the local exponents differs from that for the 
isotherm T* = 1.8: here z is a decreasing function of H, 
almost identical to the one for p = 3 (see Fig. HJb)). 



120 



100 



60 



40 



20 







l | l l l | l l l 


(a) - 
i 








~ -2 






i 


c 
1 






p = 1.5 1 

N: 


-4 




y 

i , i 


p = 2 / V 


T= 


= T 

C 


" 6 lnlHI " 4 


N, 
p\ 3 /■ // 
















■* x ^p = 50 









-0.04 



-0.03 



-0.02 
H 



-0.01 



0.75 



0.7 



0.65 



0.6 



0.55 



0.5 



T = T 

1 c 


(b) - 






\ P = l-5 








p = 50 — " 


















p= / 3 







0.01 0.02 
IHI 



0.03 



FIG. 10: (a) Wetting layer thickness £q in units of the lat- 
tice constant a as a function of the bulk magnetic field H 
calculated along the critical isotherm T* = T c * = 2.269 (path 
(4) in Fig. for L = 500, hi = 0.8, and for the ranges 
p = 1.5,2,3,4, and 50 of the boundary fields. On the right 
the curves end at the corresponding capillary condensation 
point Hc a (T, L;p). The inset shows a log-log plot of these 
data. The thin solid line indicates the slope 8/15. (b) The 
local exponents z for £o(H) (see Eq. (f20)) ). For if->0, with- 
out the occurrence of capillary condensation, i.e., for L — ► oo, 
the expected exponent is v/A = 8/15 (A). The dotted lines 
indicate the expected extrapolation for z(p > 2, H — > 0) with 
L — > oo. For p = 1.5 the behaviors of lo(H) and z(H — ► 0) 
are distinctly different from those for p > 2. 

For all decay exponents p > 2 studied for T* = 2.25 
the values of the local exponents z(H — > 0) seem to at- 
tain values between 0.48 and 0.52. They would have to 
bend upwards for H — > in order to reach the expected 
value of 8/15 ~ 0.533. For the size L = 500 capillary 
condensation does not allow to reach the ultimate com- 
plete wetting regime below the crossover line (see 



Fig.rjfor path (3)). 

For p = 3 and 2, the variation of the local exponents 
z along the critical isotherm T = T c is very similar to 
the one in Fig. [DJb). In both cases, z first decreases from 
the value rts 0.55 above the expected one 8/15 ~ 0.53 to 
the value « 0.52 at |if | « 0.001 and only then, for even 
smaller values of H, they seem to increase towards the 
expected value 8/15. In the limit H — > 0, for p = 50 
and 4 the local exponents seem to approach the expected 
value 8/15 from above whereas for p = 3 and 2 they will 
have to bend upwards do so from below. 
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FIG. 11: Adsorption F(v T J") (a) and the corresponding local 
exponents z (b) along the isotherm T* = 1.8 (path (1) in 
Fig. |3j calculated within the DMRG approach for a d — 2 
Ising strip of width L — 500 subject to symmetric boundary 
fields of strength hi — 0.8 and various ranges p corresponding 
to power law decays ~ j~ p . For small \H\ all curves end 
at H = H ca (T, L;p) where capillary condensation sets in. 
Dotted lines indicate the expected extrapolations to H = 
for semi-infinite systems: I3 c a ° = 1/3 (*) for p > 3, f3 c a ° = 1/2 
(■) for p = 2, and /3 c a ° = 2/3 (•) for p = 1.5. 
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3. Adsorption 

For the present microscopic model in the strip geome- 
try the adsorption T is defined as 

L/2 

T = Y^(m j -m b (T,H)), (23) 
i=i 

where nij = (aj) is the magnetization in row j and 
TOf,(T, H) is the corresponding bulk magnetization. In 
the limit L — > oo this definition reduces to the semi- 
infinite quantity T defined in Eq. ([3]). Along the vari- 
ous paths indicated in Fig. ([3]), at each point we have 
first calculated the bulk magnetization m& (T, H ) using 
an equivalent system but with a vanishing surface field 
in order to minimize the influence of the boundary. 
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FIG. 12: The same as in Fig. ITU along the isotherm T* = 2.25 
(path (3) in Fig. E}. In (b) {v - 0)/A = 7/15 is indicated 
(♦). The crossover to the ultimate complete wetting behavior 
is masked by the occurrence of capillary condensation. For 
p — 1.5 the behaviors of T(H) and z(H — > 0) are distinctly 
different from those for p > 2. 

Moreover, for the latter calculations as a representa- 
tive spin the one in the middle of the strip was taken, 



as it is the one least affected by finite-size effects. In or- 
der to obtain the bulk magnetization mf,(T,H) we have 
extrapolated this midpoint value m£/ 2 (T, — 0) to 

L — ► oo from data calculated for L = 300,400,500,600, 
and 700. 

The adsorption T(H) and its local exponents calcu- 
lated for a strip of width L = 500 with hi = 0.8 and 
for various ranges p of the boundary fields are shown in 
Figs. ED] [U and [13] for the paths (1), (3), and (4), re- 
spectively, as indicated in Fig. [3] 
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FIG. 13: The same as in Fig. 1121 along the critical isotherm 
T = T c (path (4) in Fig. GjJ. In (b) (y - f3)/A = 7/15 is in- 
dicated (♦). The crossover to the ultimate complete wetting 
behavior is masked by the occurrence of capillary condensa- 
tion. For p — 1.5 the behaviors of T(H) and z(H — > 0) are 
distinctly different from those for p > 2. 

According to the predictions summarized in Sec. [ill 
upon approaching the bulk coexistence line H = along 
path (1) (T* = 1.8) the adsorption is expected to diverge 
for p > 3 with the exponent (i c s ° = 1/3. As can be 
seen in Fig. ITlTb) the local exponents calculated for p = 
50, 4, and 3 tend towards this predicted value for H — > 
but do not reach it due to the occurrence of capillary 
condensation. For p — 2 the local exponents as a function 
of the bulk field exhibit a negative curvature, consistent 
with the behavior of z for the thickness £q of the wetting 
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layer (see previous subsection), and seemingly tend to 
the exponent 1/2, which agrees with the predicted non- 
universal behavior for p < 3, i.e., fi c ° = l/p. In order to 
reach the expected exponent (3 C S ° = 2/3 for p = 1.5, the 
local exponents have to vary in a non-monotonic way for 
H -> (see Fig-Htb)). 

As expected from the analysis of the thickness of the 
wetting layer, the local exponents for the adsorption 
T(H) along path (2) (T* = 2.0) behave in a similar way 
as along path (1) (T* = 1.8), i.e., for p = 50,4, and 3 
they tend to the complete wetting exponent fif = 1/3, 
whereas for p = 2 to the non-universal value /3°° — l/p — 
1/2 (data not shown). As before, the asymptotic behav- 
ior is not yet completely reached when capillary conden- 
sation occurs. We note that the theoretical considera- 
tions in Sec. [TT] do not predict any specific power law 
behavior in terms of H for the accessible H values along 
this isotherm (2). 

Paths (3) and (4) lie in the critical adsorption regime 
and one expects (see Eq. ©) T(H -► 0) ~ \H\-^~^' A 
as H — ► for p > 15/8. One can see from Figs. [T2Tb) 
and [T3T b) that for both p = 50 and p = 4 along both 
isotherms the local exponents seemingly tend to the pre- 
dicted value iy — 0)/A = 7/15 from above. In order to 
reach this value for p = 3 and 2, the local exponents have 
to bend upwards for H — > 0; whether for path (3) they 
ultimately tend to the complete wetting value (3 C S ° = 1/3 
and 1/2, respectively, after crossing both crossover lines 

Hcl and He? is unclear because for the finite values of L 
considered here the asymptotic behavior is undercut by 
capillary condensation. (Note that for f3 c ° — 1/2 one has 

= d 1 \T\ 2v .) According to Fig.[3]for T* = 2.25 this 
complete wetting regime is expected to be very narrow. 

For p = 1.5 (< 15/8) the behavior of the critical ad- 
sorption is expected to be non-universal. Indeed, as com- 
pared to the other values of p we observe a distinctly dif- 
ferent behavior of the local exponents along both paths 
(3) and (4) with z(H — > 0) strongly increasing. (Note 
that for p = 1.5 one has H& = d^T^' 2 .) 

For r ^ the capillary condensation line H ca (r, L) 
is expected (4f| to approach the bulk coexistence line 
H = as 1/L for L — > oo in accordance with the Kelvin 
equation (Eq. (fl"4|) ). Near T c , i.e., for |r| <C 1 scaling 
arguments lead to the scaling behavior (see Eqs. © and 

my- 

H ca ( Tl L) = \T\ A g{L/iM\r\-^), (24) 

where Ai is the surface gap exponent(Ai [d = 2) = 1/2) 
due to the presence of surface fields [3]; note that g < 0. 
(In Eq. ([24]) the scaling of hi holds if p is sufficiently large, 
i.e., p > 15/8.) Since £ depends on An (see Eq. ©), in 
Eq. (|2~4"|) the first scaling variable of g depends also on 
H ca (T,L), rendering an implicit equation for iJ ca (r, L). 
To leading order H ca (T 7 L) is given by Eq. ([24)1 with £ 
from Eq. (JTUJ) . In the following hi is considered to be 
large enough so that g{x,u — > oo) = g(x) < becomes 
independent of the second scaling variable. The analytic 



property g(x — > oo) ~ 1/x renders H ca (r,L — > oo) ~ 
\t\ a ~ v /L for t 7^ fixed and L — > oo, in accordance 
with the Kelvin equation. On the other hand, for L fixed 
and t — ► one has g(x — * 0) ~ x~ A ^ (see below), which 
ensures that the limit r — > renders a nontrivial function 
of L: H ca (r -> 0, L) - L~ A I V with A/v = 15/8 in d = 
2. For weak surface fields H ca (r,L) - h x \r\ A - Al g{L / '£) 
with g(x 0) ~ x-( A - Al ^ u : g < 0, so that H ca (r -> 
0,L)-i~( A " Al )/ y 0. 

The capillary condensation line ends at a capillary crit- 
ical point, at which the free energy of the confined system 
expressed in terms of the scaling variables (L/£, if |r|~ A ) 
is singular at a point (xo,yo) implying Eq. (|15p . which 
is consistent with the scaling behavior H ca (T, L — > oo) ~ 
\t\ v /L. This scaling behavior of the capillary conden- 
sation line is sketched in Fig. [21 
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FIG. 14: The log-log plot of the scaling function g(x = L/f) 
(see Eq. of the capillary condensation line H ca (r, L) cal- 

culated for several large values of L and for a short-ranged 
(p = 50) and a long-ranged (p = 3) boundary field. In the 
limit x = L\t\ v — > the scaling function saturates at a con- 
stant nonzero value whereas g(x — > oo) ~ x w with w — 0.75. 
The expected value (A — v)jv — 7/8 = 0.875 is not yet 
reached. 

With £ ~ |7"| _ly (see Eq. ©) the scaling form for the 
capillary condensation line can be written equivalently as 

H ca (r,L) = L- A ^g(L/0, (25) 

where g(x) = x A ' v g(x) < 0. We test to which extent 
scaling holds for the present model with p = 50 and 
p = 3 by plotting \H ca \L A l v versus L\r\ v for several strip 
widths L. Although we neglect the H dependence of £ 
(see above and compare Eq. ([9])) the data collapse is very 
good. From the log- log plot (see Fig.[TJJ one can see that 
g[x — ► 0) = const which implies g(x — ► 0) ~ x~ A l v (see 
above) and H ca (0,L) ~ L~ A / V . In the limit x 3> 1 our 
data fit very well to a power law g(x — * oo) ~ x w with 
the effective exponent w ~ 0.75. w is smaller than the 
value (A — v)/v = 7/8 = 0.875 which would lead to 
the predicted behavior H ca (r, L — > oo) ~ \t\ a ~ v /L. One 
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might be inclined to put the blame for the fact that w has 
not yet reached this expected asymptotic value on correc- 
tions caused by the aforementioned H dependence of the 
bulk correlation length, according to which Eq. (f2"5)l reads 
H ca {r,L)L A / v = .9( J L|rr(S(i? ca |r|- A )- 1 ) (see Eq. ©); 
this additional dependence actually spoils the scaling in 
terms of the scaling variable x = L\t\ v . Since, however, 
scaling is observed, the latter dependence must be weak 
and the too small value of the exponent w must be due 
to other corrections. 
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FIG. 15: Magnetization profiles near one wall, relative to 
their bulk values m,b(T,H), along the thermodynamic path 
(0) (see Fig. [3j shifted by a constant — 10 -6 away from the 
pseudo- phase coexistence line H ca (T, L), calculated within 
the DMRG method for d — 2 Ising strips of width L = 700 
for short-ranged ((a) p = 50) and of width L — 340 for long- 
ranged ((b) p — 2) boundary fields of strength hi — 0.8. Due 
to the slow decay, the profiles for p — 2 are presented relative 
to their values m(L/2) ^ rrib{T,H) at the midpoint of the 
film. 

Due to the scaling behavior H ca (T,L) ~ \t\ a ~ v /L, for 
L > 1 but fixed, the pseudo-capillary condensation line 
appears to approach the bulk critical point (but without 
reaching it actually due to L < oo) as |t| A-i/ = |r| 7 / 8 

which is a weaker power law than \Hcl | ~ |r| l/ / ,3s = |r| 3 



for the crossover line to complete wetting. This explains 
the observation that in d = 2, upon approaching bulk 
coexistence sufficiently close to the critical point, the 
asymptotic complete wetting regime is always preempted 
by capillary condensation. More generally, we expect 
that in confined systems the asymptotic complete wet- 
ting divergence of the adsorption can be observed only if 
A-v = 2-a-P-v> vjfif. This condition is fulfilled 
neither in d = 2 with A — v = 0.875 and v/ (3f = 3 nor 
in d = 3 with A - v = 0.936 and v j (if = iv = 1.89. 
Thus we conclude that along an isotherm close to T c the 
ultimate crossover to the complete wetting behavior can 
only be observed in systems with macroscopically large 
transverse extensions. 

In the present confined system the analogue of the ther- 
modynamic path parallel to the bulk coexistence line 
on the gas side with a small undersaturation A/i = 
const > and with T — > T c is the thermodynamic path 
-ff(o) (T, L) = H ca (T, L) + SH along the pseudo-phase co- 
existence (capillary condensation) line H ca (T,L) shifted 
slightly by SH to the spin down (a = —1) side (see path 
(0) in Fig. [3]). In view of the large computational effort 
required for the determination of the curve H ca (T, L) and 
of the value of the bulk magnetization mb(T, H) at each 
point along this line we have calculated the adsorption 
along this path only for a short-ranged boundary field 
(p = 50), for the marginal case p = 3 of long-ranged 
surface fields, and for p — 2. We have considered three 
constant shifts away from the pseudo-phase coexistence 
line: SH = -10~ 5 , -10" 6 , and -10~ 7 . 

4- Magnetization profiles 

In Fig. [13] we have plotted a selection of magnetiza- 
tion profiles m(z) along the thermodynamic path shifted 
by a constant — 10 -6 away from the pseudo-coexistence 
line H ca (T, L) for p = 50 (Fig. H5I») and for p = 2 
(Fig. [T5T b)). For p = 3 the profiles look qualitatively 
the same as for p = 2; they are both characterized by a 
much broader interfacial region and larger tails than the 
ones for p = 50. The shape of the profiles changes upon 
approaching bulk coexistence H = in a way similar to 
the case along the near-critical isotherm T* = 2.25. One 
observes a broadening of the interfacial region together 
with a rather small region where the magnetization stays 
close to the value rrib(T,H) > characteristic for being 
slightly on the upper side of the bulk coexistence curve, 
i.e., H = 0+. 



B. Thermodynamic path along the pseudo — phase 
coexistence line 

1. Thickness of the wetting layer. 

For three undersaturations SH relative to the pseudo- 
phase coexistence line Fig. ITBT a) shows the thickness of 
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the wetting layers as a function of temperature calculated 
for a fixed strip width L = 340 with strength hi = 0.8 
and ranges p = 50, 3, 2 of the boundary fields. The in- 
crease of the thickness of the wetting layer along these 
paths is very similar for all considered values of the decay 
exponent p. 
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FIG. 16: (a) Wetting layer thickness to (defined as in Figs. [9] 
and I10p as a function of temperature T* calculated for L — 
340, hi = 0.8, and various values of the exponent p, char- 
acterizing the range of the boundary fields, along the path 
(0) in Fig. [3] with constant shift 8H = — 10~ 6 away from the 
pseudo-phase coexisting line H ca (T, L; p) (see main text). For 
p = 50, in addition the temperature dependence of £o for shifts 
SH = — 10~ 5 and — 10~ 7 is shown (dotted and thick solid line, 
respectively). Symbols denote predictions from the effective 
interface Hamiltonian (EIH) for p — 50 (circle) and for p — 3 
(square) (see main text), (b) The local exponents z for £o cal- 
culated from the data with constant shift SH = — 10 -6 shown 
in (a) as a function of r = (T — T c )/T c . At the dotted and 
dashed vertical line the thermodynamic paths pass through 
the crossover line H„ and HcV , respectively, shown in Fig. [3] 
(and de facto independent of SH). Due to the finite peak of 
£o at T c , below H^) , i.e., for |r| < 0.04, the concept of the 
local exponents turns out to be less useful so that the lines in 
(b) have not been continued into this region. 

At low temperatures, for p = 2 the wetting layer 
is about three times thicker than for the short-ranged 



boundary field. Except for the close vicinity of the crit- 
ical temperature, i.e., for -0.98 < T* - T* < 1.01, 
where the finite-size effects are very pronounced, results 
obtained along paths with different shifts 8H but the 
same p collapse onto a common curve. We note that 
for p = 50 this common curve splits into three dif- 
ferent lines just above the pseudocritical temperature 
T* L (L = 340, hi = 0.8) w 2.2, i.e., when the pseudo- 
coexistence line H ca (T,L) intersects the crossover line 

He? = — c?2 1 t | A ■ For each p the thickness of the wet- 
ting film attains a certain L-dcpcndent value at the bulk 
critical temperature T c , which also depends on the value 
SH of the shift of the path (0). In order to calculate 
£q(T = T c ) for large L we have performed an extrapola- 
tion scheme for the curves £q(T) to T = T c , not taking 
into account data for temperatures higher than T* = 2.24 
because those exhibit pronounced finite-size effects be- 
yond the leading L-behavior. For p = 50 we have found 
(data not shown) that the thickness of the wetting layer 
at T = T c is proportional to L, which is expected because 
£q cannot grow larger than L/2. For p > 2 it should di- 
verge for L — > oo. Above T c we observe a rapid decrease 
of £o as function of T. 

For increasing T < T c in Fig. [TST a) we observe a con- 
tinuous increase of the thickness of the wetting layer. We 
recall that, according to effective interface Hamiltonian 
studies, sufficiently far away from criticality, i.e., for £q 3> 
£, the wetting layer thickness should be independent of 
temperature along the path of constant undersaturation 
A/i, i.e., H = const with respect to bulk coexistence (see 
the last paragraph of Sec. [TTJ) . The present results arc 
calculated along H = H {0) (L, T) = H ca (T, L) + 8H with 
SH = const < such that \H\ decreases for increasing 
temperatures, which should result in thicker films. This 
is in accordance with the data in Fig. [TBl a). In order to 
check whether this increase is captured by effective inter- 
facial models we adopt a simple effective interface poten- 
tial W(£) valid for quasi short-ranged (p = 50) boundary 
fields as well as for long-ranged boundary fields with a 
decay exponent p = 3 or 2. Taking an undersaturation H 
into account, we make the ansatz [3j (in units of /cbT c / a 2 ) 

W{£) = Ar (p - 1] +B£- K + \H\Am£ (26) 

where an effective repulsive interaction ~ B, which ac- 
counts for the gain in entropy of the unbinding inter- 
face, has to be added to the leading long-ranged po- 
tential energy term ~ A > (see Subsec. Ill A|) [261 ] . 
One has B = B \k B T '/E(T)| > H, where S(T) is 
the surface stiffness of the interface, Bq = const and 
Am = 2|m^(T)|. The exponent k describing thermal 
wandering is given by k = 2{d— l)/(3 — d), which equals 
2 in d = 2 (see Subsec. MM Ell- 
in Eq. (|26|) the second (entropic) term dominates in 
d = 2 for p > 3. Thus for short-ranged forces (p> 1) 
the equilibrium wetting film thickness varies as function 
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of T along a thermodynamic path H(T) as 

to = (i?ofcBT/(I](r)|m[ 0) (T)|)) 1 /3|^ (T )|-i/3 ; p > 3. 

(27) 

For p = 3 the entropic term competes with the interac- 
tion term. This leads to the equilibrium thickness of the 
wetting film 

io = 

((A + B k B T/V{T))/\ m f\T)\y/3\H{T)\~ l l\ p = 3. 

(28) 

A is proportional to Am = 2\m < b °\T)\ so that 

A/\m b ° (T)\ only weakly depends on T. For p = 2 
the first term in Eq. (|26p dominates and the temper- 
ature dependence of the film thickness follows that of 
H(T)~H ca (T, L): 

£o = (A/(2\m[°\T)\)) 1 / 2 \H ca (T 1 L)\- 1 / 2 , p = 2. (29) 

For the d = 2 Ising model the surface stiffness is 
known exactly as S(T)/(fc s T) = sinh2(A' - if*) where 
A" = J/{k B T) and tanhA* = exp(-2if). Both E(T) 
and Ato(T) = 2|m[, 0) (r)| decrease with increasing tem- 
perature. Therefore, for all values of p, £q is expected to 
increase along the pseudo-coexistence line H ca (T,L). 

Sufficiently close to the bulk critical point the sur- 
face stiffness vanishes with the same power law as the 
surface tension of the interface a, i.e., S(r — > 0) ~ 
a ~ \ T \( d - l > = \t\ for the d = 2 Ising model. Since 
Am ~ \t\& it follows that in the critical region for short- 
ranged forces 

4Hrr ( ^ (d - 1M/3 iiMT,£)r 1/3 = 

\T\- 3 / s \H ca (T,L)\-V*, p = 3. 

(30) 

with H cci (t,L) given by Eq. ([25"f . 

Close to T c and in the temperature range for which 
the pseudo-coexistence line H ca (r, L) lies in between the 

crossover lines H^r = — di\r\ 3L/ and Hc^ = — d2|r| A we 
have found (see Fig. [14]), for both p = 50 and p = 3, that 
the absolute value of the scaling function g(x = L/tj) of 
the pseudo-coexistence line (see Eq. ([23)) ) decreases and 
crosses over from the power law behavior g(x — ► oo) ~ x w 
with w ~ 0.75 for |r| ^ to a constant value g{x = 0) 
at t = 0. This implies that, for decreasing |r|, 4)( T ) 
crosses over from an increase ~ |r| -5 / 8 to one ~ |-r | 3 / 8 
(for L fixed). Thus the effective interface model predicts 
that, along the path (0), £${t — ► 0) is an increasing func- 
tion of temperature. But in view of the aforementioned 
crossover one cannot expect a purely algebraic behavior 
in this region. 

For 12 < L\t\ v < 50, the variation of the pseudo- 
coexistence line H ca (T,L) at fixed L is approximately 
algebraic with an effective exponent w ~ 0.75. Insert- 
ing H = H ca (T,L) ~ |r| ' 75 into the expression for £q 



(Eq. d3DJ|) gives £ - |r|-°- 625 for both p = 50 and p = 3. 
For L\t\ v <C 1 one has H ca — > const, so that in this limit 
Eq. (301) predicts £ (t ->0)~ |t|- 3 / 8 . 

For p = 2 the pseudo-coexistence line H ca (T, L) seems 
to follow a power law with the same effective exponent 
as for p = 50 and 3 for reduced temperatures smaller 
than the crossover line H cr (l) (see Fig. Therefore on 
the basis of Eq. (29]) for 12 < L\t\ v < 50 one expects 
£ Q ~ |rr - 375 . 

For p = 50 the behavior of the local exponents z de- 
termined from the numerical data obtained for L = 340 
and shown in Fig. 116( a) is consistent with the above pre- 
dictions (see Fig. ITBT b)). We find that in between the 
crossover lines H^) = — c?i|r| 3ly and H„ = —d 2 \T\ A the 
local exponents z decrease continuously upon decreasing 
|r|. For 0.12 < L\t\ v < 50, which with L = 340 and 
v = 1 corresponds to 0.04 < |t| < 0.15, the value of z 
for p = 50 changes from ~ 0.80 to ~ 0.52. No satura- 
tion at the value 0.625, which would reflect the power 
law g(x) ~ x 75 (see Eq. (2"5"]) and the discussion above), 
is observed. However, the above temperature range lies 
outside the range of validity of the asymptotic power law 
behaviors of £(T) and Am(T) leading to deviations from 
the algebraic variation £ ~ | T | -0.625 cji scussec j above. 

In order to test the validity of the effective inter- 
face Hamiltonian approach we compare our full DMRG 
data with the predictions for the increase of the wet- 
ting film thickness along the line of pseudo-phase co- 
existence H ca (T 7 L) stemming from Eq. ([2"T]) by substi- 
tuting our numerical data for H(T) = H ca (T,L) and 
Am(T) = 2|m&(T)| and by using the analytic expression 
for E/(fcsT) given above. Treating Bq as a fitting pa- 
rameter we find very good agreement between these two 
approaches for Bq = 0.729 (see Fig. [TBTa)) in the wide 
temperature range below T* ~ 2.215. We note, that for 

higher temperatures the system crosses over to the criti- 

(2) 

cal adsorption regime (beyond the crossover line mP). 

For p = 3 and 2 we observe that at the local 

exponents z saturate at the values ~ 0.55 and ~ 0.44, 
respectively. In the case p = 3 the exponent is slightly 
smaller than the prediction 0.625 based on the scaling 
behavior obtained within the effective interface model, 
whereas for p = 2 it is slightly larger than the value 
0.375 expected also within the effective interface model. 

We repeat the procedure described above to check the 
agreement of the effective interface Hamiltonian predic- 
tion in Eq. (28]) with our DMRG data for p = 3. We 
adopt Bq = 0.729, i.e., the value determined from the 
fitting of Eq. (27]) to the DMRG results for p = 50, and 
treat A/\rrib(T)\ f» const = C as a free parameter. Very 
good agreement is obtained for C ~ 1.5 and tempera- 
tures T* < 2.15 (see Fig. HBTa)). 

However, the assumption A/\rrib(T)\ ~ const does not 
lead to an agreement of the DMRG data for p = 2 (not 
shown) with the curve predicted by Eq. (|29[) with H(T) = 
H ca (T, L). Thus we conclude that in the case p = 2 (not 
shown) the simple effective interface Hamiltonian given 
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by Eq. (HI fails to describe the full DMRG data. 



2. Adsorption 

Results for the adsorption obtained for L = 340, h\ = 
0.8, and p — 50, 3, 2 are shown as a function of T in 
Fig-HKa). 
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FIG. 17: (a) Adsorption T as a function of T* along the 
path (0) in Fig. [3] with a constant shift SH = — 1CP 6 away 
from the pseudo-phase coexisting line H ca (T, L;p) (see main 
text), calculated within the DMRG method for d = 2 Ising 
strips of width L = 340 for short-ranged (p = 50) and long- 
ranged (p — 3, 2) boundary fields of strength hi = 0.8. For 
p — 50 the temperature dependence of F for shifts SH = 
— 10 -5 and — 10~ 7 is also shown (dotted and thick solid line, 
respectively), (b) The local exponents z for F calculated from 
the data shown in (a) for the constant shift SH = — 10 -6 . At 
the dotted and dashed vertical line the thermodynamic path 
passes through the crossover lines H^) and H 1 ^ , respectively, 
shown in Fig. [3] (and de facto independent of SH). As in 
Fig. 116T b V in (b) the lines do not continue into the region 
\t\ < 0.04. 

In accordance with the expectations, the behavior of 
r(T) upon increasing the temperature is qualitatively 
similar to the behavior of £q(T) in Fig. [Tffla). Results 



for different shifts 5H form a common curve above T c 
and for T* < 2.2. Close to T c the adsorption exhibits a 
maximum which signals the vicinity of the end of the 
pseudo-phase coexistence line, i.e., the pseudo-critical 
point of a finite system. The position of this maximum 
shifts towards T c for wider strips. At T c , T attains a 
finite value T(r = 0, L). As expected, for long-ranged 
boundary fields the adsorption is slightly stronger than 
for short-ranged forces. 

According to the discussion in the last paragraph of 
Sec. HH for p > 3 and in between the two crossover lines 
= -di|r| 3l/ and H$ = -d 2 \r\ A (sec Fig. |3) T is 
expected to increase ~ \t\~^ u ~P\ which for d = 2 cor- 
responds to an exponent -0.875, provided the variation 
of £ is dominated by its dependence on r (Eq. I|10p). 
This is the case along paths of constant undersaturation 
An = H = const. Along the pseudo-coexistence line 
and for p = 50 the local exponents of T vary strongly 
with temperature, also below the crossover line (sec 
Fig- fTTTb)). This is not compatible with the aforemen- 
tioned algebraic variation of T for H = const. In view 
of the discussion of the corresponding behavior of the 
wetting film thickness in the previous subsection, this 
difference can be traced back to the fact that the pseudo- 
coexistence line does not run in parallel to the bulk co- 
existence line. The decrease of the absolute value of the 
bulk field H ca (T,L) upon increasing temperature mat- 
ters and in the scaling analysis of T the whole scaling 
form of £ (Eq. ©) has to be taken into account. Accord- 
ing to Eqs. and (O this is encoded in the product of 
the scaling functions H± and K± as function of -ff|r|~ A , 
which alters the temperature dependence ~ \t\~^ u ~^\ 
Due to r ~ ioAm and the fact that Am decreases upon 
decreasing the reduced temperature |r|, one expects that 
the increase of T is weaker than the increase of £q . This 
is consistent with our numerical data for the local expo- 
nents z. At each point along the path (0) the value of z 
obtained for T is smaller than the corresponding value of 
z obtained for Iq (see Figs. HBT b) and [TTTb)). 

For p = 3 and p = 2 the variation of the local exponents 
of the adsorption in the temperature range in between 
the crossover lines and He? is less pronounced than 
that for p = 50. Near Her we find saturation for the 
local exponents z at values which are smaller than the 
saturation values of z of the corresponding thickness £q of 
the wetting film, i.e., z ~ 0.49 and z ~ 0.36 for p = 3 and 
p = 2, respectively (see Fig. [TTl b)). This is consistent 
with the relation T ~ igAm and can be explained by 
the same arguments as used in the previous subsection. 
We note, that Am - \t\ , with /3 = 1/8 = 0.125 for the 
d = 2 Ising model, leads to bigger differences in the local 
exponents for £q and T than the aforementioned ones 
observed in the range of their saturation [p = 3: z(£q) — 
z(T) = 0.55 - 0.49 = 0.06 < 0.125; p = 2: z(£ ) - z{T) = 
0.44 - 0.36 = 0.08 < 0.125]. However, the above power 
laws are satisfied only sufficiently close to T c , which is not 
reached along H ca (T, L) for the finite values of L studied 
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here. 



VII. SUMMARY AND CONCLUSIONS 

We have studied two-dimensional Ising ferromagncts 
in strip geometries of width L and with long-ranged 
boundary fields (V? = ! ^,h 1 > 0) (Eqs. JT2]| and 
(fT3|) ). Based on scaling theory and the density-matrix 
renormalization-group method we have obtained the fol- 
lowing main results: 

(1) In Sects . HI and HT1 we have discussed the theoretical 
framework and the expectations due to scaling theory for 
the interplay between critical adsorption and complete 
wetting, including the crossover between them. For semi- 
infinite systems the corresponding scaling arguments lead 
to the schematic phase diagram shown in Fig. [1] 

(2) The location of the pseudo-phase coexistence line 
H ca (T, L;p) of capillary condensation has been deter- 
mined for various ranges p of the boundary fields (Fig. [3]). 
For positive and parallel surface fields capillary conden- 
sation occurs at negative values H ca (T, L;p) of the bulk 
field H . Increasing the range of the boundary fields shifts 
the pseudo-phase coexistence line towards more negative 
values of the bulk field H (see Figs. [3] and 2]). Relative to 
the pseudo-capillary condensation line in the case that 
the surface fields act only on the two boundary layers 
j = 1 and j = L, i.e., for p = oo, H ca (T, L;p) varies 
exponentially as function of the range p (Fig. [3J . 

(3) A critical wetting transition has been identified for 
p = 2, i.e., for a boundary field decaying slower than 
the marginal one (p = 3) (sec Fig. [SJ. These findings 
disprove previous claims by Kroll and Lipowsky [37j that 
in d — 2 for p < 3 there is no critical wetting transition 
at finite values of hi. 

(4) Within the accessible range of values for T and 
H the magnetization profiles both along isotherms (see 
Figs. [5] and [7]) and along the pseudo-capillary conden- 
sation line (see Fig. [3]) are not slab-like. They exhibit a 
wide interfacial region and significant tails. Primarily the 
profiles approach their plateau values in the middle of the 
strip exponentially as function of the distance from the 
wall. The width of the emerging interfacial region and 
the decay length of the tails are proportional to the bulk 
correlation length and thus grow upon approaching T c or 
the pseudo-capillary condensation line. In the presence of 
long-ranged boundary fields the exponential decay of the 
profiles towards their bulk values mi,(H, T) is followed by 
an algebraic decay j~ p which finally is distorted by the 
presence of the distant wall. Features similar to those 
of these order parameter profiles have been inferred from 
neutron rcflcctrometry for the composition profile of a 
wetting film in a binary liquid mixture of n-hexane and 
pcrfluoro-n-hexane [47| . 

(5) The variation of the thickness £q of the wetting 
layer along various isotherms (Figs. [5] fTU|) and along the 
pseudo-phase coexistence line (Fig. ll6[) has been analyzed 
for wide strips and different ranges p of the boundary 



field. Along both types of path we have found a grad- 
ual increase of £q upon approaching T c . The asymptotic 
divergence of £o along the isotherms is preempted by cap- 
illary condensation. As discussed on general grounds at 
the end of Subsec. IVI A 3[ both in d = 2 and d = 3 
along the isotherm close to T c the ultimate crossover to 
the complete wetting behavior (see hatched region and 
path II in Fig. [T|) can only be observed in systems with 
macroscopically large transverse extensions. Along the 
isotherms and within the complete wetting regime the 
suitably defined local exponents of £q (Eq. (f2"0"|) ) tend to 
approach the values predicted by the corresponding ef- 
fective interface Hamiltonian, i.e., j3f = 1/3 for p > 3 
and (3 C S ° = 1/p for p < 3 (Figs. |5] and EJ). These re- 
sults further support our findings that a critical wetting 
transition exists also for the values p = 2 and 1.5 of the 
decay exponent of the boundary field, i.e., smaller than 
the marginal case p = 3. Along the pseudo-phase co- 
existence line H ca (T,L) the variation of the thickness of 
the wetting film is determined by the functional form of 
the temperature dependence of H ca (T, L) and for p > 2 
it agrees with the predictions of an effective interface 
Hamiltonian (Eq. (j2~7| and Fig. [TBf a)). For long-ranged 
boundary fields with p — 3 and p = 2, the increase of £q 
within a certain range of temperature can be described 
by a power law with an effective exponent which is, how- 
ever, not universal. It varies from ca 0.55 for p = 3 to ca 
0.44 for p = 2 (Fig.QStb)). 

(6) For p = 50, 4, 3, and 2 the adsorption V (Eqs. © 
and l|23p) has been calculated along various isotherms 
(see Figs. fTT1fT5|) . Along these thermodynamic paths, 
over a wide range of the bulk field H the adsorption 
r exhibits a continuous increase upon approaching the 
pseudo-phase coexistence line H ca (T, L;p). Within the 
accessible range of values for H this increase cannot 
be described by simple power laws. Only very close 
to H ca (T, L;p) the local exponents of T start to ap- 
proach their predicted values, i.e., (3f(p > 3) = 1/3 and 
fts°(P = 2) = 1/2 for complete wetting (see Subsec. Ill Al 
and Figs. [H and Q2)) and (u - /3)/A = 7/15 for crit- 
ical adsorption (see Eq. Subsec. Ill Bt and Fig. [T3|) . 
However, these results are not entirely conclusive because 
the asymptotic regimes are preempted by capillary con- 
densation so that one is unable to detect the crossover 
from critical adsorption to complete wetting which is ex- 
pected to occur sufficiently close to bulk coexistence for 
isotherms with T close to T c . For the low temperature 
isotherms, which lie entirely within the complete wetting 
regime, the cases p = 1.5,2, i.e., for p smaller than the 
marginal value p = 3, differ distinctly from the cases 
p > 3 (sec Subsec. Ill Al and Fig. fTTj) . which is a clear in- 
dication of the non-universality of complete wetting for 
these cases. For p — 1.5 the long-ranged part of the 
boundary field is relevant in the RG sense and the be- 
havior of r is non-universal along the critical isotherm 
(see Fig.H3Ib)). 

(7) The finite-size scaling predictions for the capillary 
condensation line (see Eqs. (|2"4"|) and (f2"S"|) as well as Figs. [2] 
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and [14]) arc satisfied for both short-ranged (p = 50) and 
long-ranged (p = 3) boundary fields. 

(8) For p = 50, 3, and 2 the magnetization pro- 
files (Fig. [T5"|) and the adsorption T (Fig. ITT)) have been 
calculated along the pseudo-capillary condensation line 
H ca (T, L). Over a wide range of temperatures below T c , 
r increases gradually (sec Fig. [T7|) . Near criticality its in- 
crease cannot, however, be described by a simple power 
law behavior. The reason is that along this thermody- 
namic path the temperature dependence of T is deter- 
mined also by the temperature variation of H ca (T, L), 
similarly to the thickness £q of the wetting film which is 
approximately related to T as T ~ io^m where Am(T) is 
the temperature dependent difference between the mag- 
netization of the spin up and the spin down bulk phases. 
Our theoretical results cannot be directly compared with 
the experimental findings in Rcfs. 1 151 and [16[, in which 
the thickness of the wetting layer [15| and the adsorption 
[l(| were measured along the paths corresponding to the 
path (1) in Fig.[TJ i.e., in these experiments the capillary 
condensation was absent. The data obtained from both 
experiments show that the behavior of Iq and T can be 
described by a power law but with an effective exponent 
which differs from both that for complete wetting and 
that for critical adsorption. Similarly, in our theoretical 



analysis of two-dimensional Ising strips, for long-ranged 
boundary fields with p = 3 and p = 2 the increase of F, 
for still sufficiently large reduced temperatures |r|, can 
be described by a power law with an effective exponent, 
here with a value of ca. 0.44 for p = 3 and of ca. 0.36 for 
p = 2 (Fig. ED. 

(9) Our results indicate that the asymptotic behavior 
for H — > occurs in a much narrower regime than one 
would expect and therefore requires much larger system 
sizes to suppress capillary condensation sufficiently in or- 
der to reach the asymptotic regime. One indication for 
having entered the asymptotic regime could be that a 
fully developed wetting layer has been formed. For the 
system sizes L and thus the values of H accessible in our 
calculations this does not yet fully occur. 
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